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Abstract — In OFDM systems, the symbol timing synchronization 
is needed to find an estimate of where the symbol starts. The most 
popular of the pilot-aided timing estimators is the one proposed 
by Schmidl and Cox (S&C) ["Robust Frequency and Timing Syn- 
chronization for OFDM," IEEE Transactions on Communications, 
vol. 45, no. 12, pp. 1613-1621, Dec. 1997]. Generally, the proper 
action of the timing estimator may be strongly affected by the pres- 
ence of disturbances, like narrowband interference (NBI) signals. 
This paper investigates the effect of NBI signals, which may arise 
in the OFDM band as the spectrum becomes more crowded, on the 
timing synchronization for OFDM systems. The performance of 
the S&C symbol timing synchronizer is evaluated in an analytical 
way in the presence of narrowband interferers. Also, the probabil- 
ities of both missing and false detection of a training sequence are 
addressed. Further, simulations have been carried out to verity the 
validity of the approximations in the analysis. 

Index Terms — Interference, orthogonal-f^quency-division 
multiplexing (OFDM), symbol timing synchronization, training 
symbol. 



I. Introduction 

ONE of the most challenging issues for fumre wireless 
communications is the provision of new spectral re- 
sources and the more efficient use of existing frequency bands. 
Recently, new spectrum management strategies have been 
proposed [2], [3] which allow two or more systems to share 
the same frequency band. The broadband very high frequency 
(B-VHF) project! [4], [5] which aims to develop a new inte- 
grated broadband VHF system for aeronautical voice and data 
hnk communications based on multicarrier (MC) technology, 
is a good example of an overlay system. In this project, the 
MC system is intended to share the parts of the VHF spectrum 
with currenfly existing narrowband (NB) systems such as voice 
double side band/amplitude modulation (DSB-AM) signals 
and VHF digital links. These NB systems are considered as 
interference to the MC system and limit the effectiveness of 
the multicarrier system. According to the authors' knowledge, 
no previous research was done to investigate the effect of 
these NB signals on the performance of a MC system. In [6], 
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Coulson investigated the effect of narrowband interference 
(NBI) on pilot symbol detection and frequency synchronization 
for OFDM systems; tiie NBI model in [6] was an unmodulated 
complex sinusoid. Alfliough this model is an accurate model for 
different interference sources [7] such as carrier feed-tiirough, 
baby monitors, and garage doors, this model is not sufficient 
to describe interference sources with a bandwidth that is larger 
than a small fraction of the carrier spacing of the MC system. 
This motivates us to study the effect of digitally modulated NBI 
signals on the OFDM systems. In [8], we have investigated 
the effect of digitally modulated NBI signals on frequency 
ambiguity resolution of OFDM systems. The result of [8] indi- 
cates that the performance of the integer frequency estimator is 
strongly affected by NBI. In this paper, we evaluate die impact 
of digitally modulated NBI signals on timing synchronization, 
which is one of the most important parts of the OFDM receiver. 
The proposed interference model for NBI can also be used to 
describe the interference in other applications where a MC 
system has to coexist with NB signals. 

In symbol timing synchronization, the main objective is to 
know when a received OFDM symbol starts. In the case of an 
ideal channel, the synchronizer must select a timing widiin the 
range of the cyclic prefix to avoid intercarrier interference (ICI) 
and intersymbol interference (ISI) [9]-[ll]. Indeed, a timing 
offset not exceeding the length of the cyclic prefix gives only 
rise to a carrier-dependent phase rotation [10], which can easily 
be corrected by die channel equalizer at die FFT outputs. How- 
ever, a time shift larger tiian the cyclic prefix will cause die in- 
noduction of ICI and ISI. In a dispersive channel, the interval 
over which die timing may range is reduced as compared to the 
case of the ideal channel, with an amount equal to the impulse 
response duration of the channel. 

Symbol timing synchronization at the receiver can be 
achieved by inserting some synchronization symbols at the 
transmitter as pilot OFDM symbols [1], [12], [13]. These 
symbols are then picked up by the receiver to synchronize the 
symbol clock. The most popular of the pilot-aided algorithms 
is the method proposed by Schmidl and Cox (S&C) [1]. The 
method relies on searching for a pilot OFDM symbol with two 
identical halves in time domain. This symbol is generated by 
transmitting pseudonoise (PN) data symbols at even frequen- 
cies, while zeros are used on the odd frequencies. The S&C 
estimator provides a simple and robust estimation for symbol 
timing. 

In this paper, we investigate the performance of the S&C 
method in the presence of NBI signals. The paper is organized 
as follows. In Section II, we describe the OFDM system in the 
presence of NBI signals. In Section HI, the statistical properties 
of the timing metric of the S&C estimator in the presence of 
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Fig. 1. Block diagram of OFDM system including interfering signals. 



interference are derived. The simulation and analytical results 
are discussed in Section IV. Finally, conclusions are given in 
Section V. 

II. System Description 

The basic block diagram of the OFDM system and NBI sig- 
nals is shown in Fig. 1 . In the OFDM transmitter, the data stream 
is grouped in blocks of Nu data symbols. Next, an N- point in- 
verse fast Fourier transform (IFFT) is performed on each data 
block, where AT > iV„ and a cyclic prefix of length v is inserted. 
The nth time-domain sample of the ith data carrying OFDM 
block can be written as 

- ^ 4 e^'^'^*"/^ -v<n<N-l (1) 




Fig. 2. Spectrum of baseband OFDM and one interfering signal. 



where /„ is the set of Nu carrier indexes and a\ is the fcth data 
symbol of the «th OFDM block; data symbols are independently 
and identically distributed (i.i.d.) with zero mean and variance 
[hlH = ^s- The M data OFDM blocks are preceded by 
one OFDM training symbol; the training symbol consists of two 
identical halves in the time domain. The time-domain samples 
of the training symbol can be written as in (1), but with the 
symbols modulated on the odd carriers equal to zero and on even 
carriers equal to a pseudonoise sequence; the energy per symbol 
on the even carriers equals i2Es). The time-domain signal of the 
baseband OFDM signal Su{t) consists of the concatenation of 
the two pilot block and the M data OFDM blocks, as follows: 

M N-1 

Suit) = E E ''i'')Po(t - - i{N + (2) 

where |i = 0 corresponds to the pilot block and i = 1, . . . , M 
corresponds to the M data carrying blocks, po(t) is the transmit 
pulse of the OFDM system, and I/Tq is the sample rate. The 
transmit pulse po (i) is assumed to be a square-root-raised cosine 
filter with roll-off qq. The baseband signal (2) is upconverted to 
the radio frequency /q. At the receiver, the signal is first down- 
converted, then fed to the matched filter and finally sampled at 
rate I/Tq, resulting in the samples ruimTo). Note that, when 
the number of carriers N is large, the sample s'(n) consists of 
a large number of contributions. Hence, taking into account the 



central limit theorem, the real and imaginary parts of s'(n) can 
be modeled as Gaussian random variables with zero mean and 
variance <jI equal to [E^ • Nu/2)/{N + u). The OFDM signal 
is disturbed by additive white Gaussian noise (AWGN) with un- 
correlated real and imaginary parts, each having variance a"^- 
The signal-to-noise ratio (SNR) at the output of the matched 
filter is defined as a^/a^. 

Further, the signal is disturbed by narrowband interference re- 
siding within the same frequency band as the wideband OFDM 
signal as shown in Fig. 2. The interfering signal may be 
modeled as the sum of Nj narrowband interfering signals 



where sj(t) is the hh NBI component 



(3) 



si{t)= ^ bi,hPi{t-hTi-Ti)-e^^<f''+f-'^' (4) 

where 6;,^ is the hth interfering data symbol of the lih interferer, 
pi{t) is the transmit pulse of the Ith interferer, r; is its delay, and 
1/Ti its sample rate; pi{t) is assumed to be a square-root-raised 
cosine filter with roll-off ai. The interfering signal is modulated 
to radio frequency fo + fc,u where /c,j is the carrier frequency 
deviation from /o for the Zth interferer. The total NBI signal may 
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at the output of the matched filter of the OFDM receiver 



^/W = E E hh^'''''^"'''''9i{t-hTi) (5) 

where gi{t) is the convolution of po(*) and pi{t — 77) 
exp(j27r/c,if). It is assumed that the interfering symbols are 
I with each other, i.e., E^i^hbii^h'] = E{SwShh', 
where liS; is the energy fJer symbol of the Zth interferer. Further, 
the intjrfering data symbols bi^h are statistically independent 
of the (])FDM data symbols aj.. The signal-to-interference ratio 
(SIR) £ t the input of the receiver is defined as 
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Statistical Properties of the Timing Metric 

IN THE Presence of NBI 
e S&C algorithm [1], the training symbol contains two 
halves in the time domain and the timing delay esti- 
earches for the peak of correlation between the matched 
utput samples that are separated by half an OFDM 
The sample r(mTo) at the output of the matched filter 
OFDM receiver consists of a useful signal r„(mro), 
component nimTo), and a noise component 
It may be expressed as 
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present, are as small as possible. To find these probabilities, 
and hence the threshold, we need to know the statistical prop- 
erties of the timing metric. In [1], the statistical properties of 
the timing metric at the optimum timing point and a position 
outside the training sequence were investigated. In this section, 
we extend the results of [1] to obtain the statistical properties 
of the timing metric in the presence of NBI. 

First, we determine the mean of the timing metric in 
the presence of NBI signals as function of the timing po- 
sition d. For sufficiently large values of SIR and SNR, 
E [\P{d)f] ^ \E [P(d)]|^ V - (iV/2 + jy) < d < N/2-1: 
the imaginary part of P(d) is small as compared with its real 
part and will be neglected. As the variances of |P(rf)p and 
R{d)^ are small, E[M{d)] w \E[P{d)]\^/E [{Rid)f]. so 
that 

E[M{d)] = l\!l:^^L, -u<d<Q 

(11) 

where t/i and G are given by 



r(mTo) = uimTo) + ri{mTo) + w(mTo). 



mbol timing estimator takes the instant d, where the 
metric 



M{d) 



- R{d)^ 



, as the starting point of the frame. In (8), P{d) and 
•e given by 

N/2-1 N\ \ 

= Y1 r*i{d + m)To)-r(^[d + m+jjToj (9) 



This timing metric is not only used to find the optimum timing 
instant: but also to determine whether or not a training sequence 
is rece ved. To do this, we use a threshold: when the timing 
metric exceeds this threshold for at least one value of d, we 
decide !that it is possible to detect a training sequence, whereas 
when the threshold is never exceeded, we decide that it is 
not possible to detect a training sequence. The value of the 
threshold must be selected such that the probability of missing 
a training symbol when there is one present, and the probability 
of falsely detecting a training sequence when there is none 



and i4('"''"')(x, y) is defined as 



2' 2 J 



(12) 
(13) 



(7) 



A('"-™')(x,3/) = E E Ei-gt{{m + d + x)To-hTi) 

1=1 h=-oo 

•gii{m,' + d + y)To-hTi). (14) 

When d is outside the region considered in (11), i.e., when 
d>N/2oTd<-iy- N/2, the app«)ximation E [\P{d)f] 
\E [Pid)]]^ is no longer valid. However, it can easily be shown 
that in that case E [M(d)] ~ D^/G^ (see also Section ffl-B) , 
where is given by 

= f {2al+2a?,)'+"^\^-'"H0, 0) • , f) • 

(15) 

In Section IV, we will evaluate the validity of the used ap- 
proximations by means of simulations. 

A. Optimum Timing Point 

Let us define q = \Pidopt)\/ R{dopt) as the square root of 
the timing metric M((/opt) at the optimum timing point dopt- 
To simplify the notation, we drop the subscript "opt." Due 
to the synnmetry of the pilot symbol, the signal component 
Emi'o"' Kiim + d)To) ■ r„((m + d + N/2)To) of P{d) is 
real valued. It can easily be shown that the other components 
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of P{d) have zero mean and small variance when the SIR and 
SNR are sufficiently large. Therefore, at high values of SNR 
and SIR, the imaginary part of P{d) can be neglected. From 
(7), (9), and (10), it can easily be verified that \P{d)\ and R{d) 
contain a common term a. Hence, we can rewrite \P{d)\ and 
R{d) as \P{d)\ = a + /3andi?(d) = a + 7, where a, /3, and 7, 
as well as their statistical properties, are given in the Appendix. 

As can be observed in the Appendix, the expressions for a, 
13, and 7 consist of a large number of contributions when N 
is large. Therefore, according to the central limit theorem, the 
variables a, /3, and 7 have approximate Gaussian distributions. 
Taking into account that for both the numerator and denominator 
of q, standard deviations are much smaller than the averages, q 
can be approximated by a Gaussian variable [14] with mean and 
variance given by the following approximations: 



M9« 



(16) 
(17) 



where ^i^ and <t^ is the mean and variance of x, respectively. 
In (17), we have used the approximation a(l + 6)/(l + c) cs 
o( 1 + 6 - c) as in [ 1 ] . Taking into account that q is approximately 
a Gaussian random variable, and M(d) is the square of q, it 
follows that M(d) also is approximately a Gaussian variable 
[14]: 

M{d) = (/xe, + iV (0 , <72))2 « ^2 + 2 . • iV (0 , al) (18) 

where N{x, y) is a Gaussian random variable with mean x and 
variance y. The mean of M(ti), i.e., /i^m. can be easily computed 
calculated by using (16), as follows: 

/ 1 ^' V 



where Tj and 4>i are given as 



T(=Re ^ E 9i{{m + d)To-hTi) 



N/2-1 00 I , 



(20) 
(21) 



Note that (19) corresponds to the result in (1 1) for d e 0 ] . 
Fig. 3 shows T; and 0; for one interferer {Nj — 1) as 
function of normalized interference frequency deviation, 
r = + ao)/2ro). From Fig. 3, it follows that 

T( 0;. It can easily be shown that when varying the normal- 
ized interference bandwidth (NBW) = Bi/Bq, where Bi is the 
bandwidth of l\h interferer and Bo is the bandwidth of OFDM 
signal; the same conclusions yield. Hence, the effect of the 




Normalized interference frequency, r 
Fig. 3. T, and (NBW = 0.0244, AT = 1024, ao = .25. a, = .5). 

interference on the numerator of (19) will be much smaller than 
its effect on the denominator. Therefore, decreasing the SIR 
will result in a decrease of fiM- When the interferer is located 
completely in the flat region of OFDM spectrum (see Fig. 2), 
i.e., |/c,j + (1 + a,)/2Ti\ < ((1 - ao)/2To), 9i{t) equals 

(22) 



where C is constant. Using (22) in (20), T; may be rewritten as 

{iV/2-1 00 
ei2./.,(iV/2)To ^ £ 

p^{{m + d)To-hT,-Ti) 



9i{t) = C-piit-T,)-e^^f-'* 



■PI 



(23) 

Because of the factor ei27r/..,(N/2)T„ i„ (23), it follows that Tj 
is a periodic function of fc,i when the interferer is located in the 
flat region of the OFDM spectrum. Similarly, when we substi- 
tute (22) in (2 1 ), it turns out that becomes independent of fc,i 
in the flat region of the OFDM spectrum. This can be observed 
in Fig. 3. 

According to (18), the variance of M{d) is given by 



(24) 



where can be easily computed by using the (16) and (17). 
Note that, when SIR tends to infinity, the same expression for 
is obtained as in [1]. 

B. Timing Position Outside of Training Sequence 

In this section, the statistical properties of a timing position 
d outside of the training sequence are considered when NBI is 
present. Using a similar analysis as in [1], it can easily be ver- 
ified that for sufficiently large SNR and SIR, the timing metric 
is approximated by a chi-square-distributed random variable, as 
follows: 

D2 



M(doutside) = 



2G2 



(25) 
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Fig. 4. Probability distributions of timing metric. 



where x^(2) is a chi-square-distributed random variable with 
two degrees of freedom and its mean and variance equal to 2 and 
4, respectively; G and 2? are given in (13) and (15), respectively. 
Therefore, the mean and the variance of M(rfoutside) are given 
by 

MiW(d<.„uid,) = ^ (26) 



Note that these results are the same as in [1] when SIR tends to 
infinity. 

C. Probability of Missing/False Detection of 
Training Sequence 

In Fig. 4, the probability distributions of the timing metric are 
shown at the optimum timing instant and at an instant outside 
the training symbol as calculated in Sections III- A and III-B. 
These distributions can be used to determine both the proba- 
biUty of not detecting a training sequence when one is present, 
i.e., the probability of miss (F„,iss), and the probability Pfaise of 
falsely detecting a training symbol when there is none present at 
a certain threshold value (A). From the statistical communica- 
tion theory, Pmiss(A) and the probability PfaUe(A) are calculated 
as follows: 

where erfc(a;)=2/(v'7r) /~ e'^^dz. Accordingly, Pfaise(A) = 
1 - P(A), where P(A) is the cumulative distribution function 
of a chi-square-distributed random variable with two degrees of 
freedom. The closed form of P(A) is obtained as [15] 

P(A) = 1 - e-^/2<r= (29) 

where is the variance of i.i.d. Gaussian random variables, 
which equal D^/2G^ as in (25). Hence, Pfaise(A) can be written 
as 

Pfalse(A) = e-^«VX^\ (30) 




-1 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1 

Timing offset c/( symbols) 



Fig. 5. Average timing metric at SNR = 10 dB, NBW = 0.0244, N, = 1, 
and r = 0. 



IV, Analytical and Simulation Results 
The numerical and simulation results in this paper are ob- 
tained with the following OFDM and interference parameters.^ 

• The total number of subcarriers is iV = 1024. 

• The total number of active subcarriers is iV„ = 1000. 

• The guard interval is set to about 10% of OFDM block size, 
u = 102. 

• The bandwidth of the OFDM spectrum equals 
Bo = 1024 kHz. 

• We use QPSK modulation for the data symbols of the 
OFDM and the interferer signals. 

• Transmit filters are square-root raised-cosine filters with 
roll-off factors ao = 0.25 and a; = 0.5 for OFDM and 
interfering signals, respectively. 

• Time delay of the interferers rj = 0. 

In [1], it is shown that the S&C timing algorithm is robust 
to fading channels. Therefore, we restrict our attention to an 
AWGN channel for simulation, to clearly show the effect of NBI 
on the performance of the synchronizer. However, the results of 
this paper can straightforwardly be extended to other types of 
channels. 

In Fig. 5, the average timing metric based on the obtained 
analytical expression (11) is shown as function of the timing 
oEfset.3 As can be observed, the timing metric shows a plateau 
at the maximimi. This is explained as the cyclic prefix extends 
the symmetry of the OFDM pilot synchronization symbol such 
that all timing instants within the cyclic prefix will give rise 
to the same timing metric. This plateau will result in an ambi- 
guity in the selection of the optimum timing. However, a timing 
offset within the cyclic prefix will cause no interference but only 

^For the sake of comparison between our results and results presented in [1], 
we use tlie same parameters as in [1]. However, in the actual hardware imple- 
mentation, the number of null subcarriers would be higher to ease filtering re- 
quirements. Moreover, the analysis of the paper is valid for other antialiasing 
filters than square-root raised-cosine filter; one just has to plug in the correct 
expression for these filters in our equations. 

Ht is easy to verify that if the normalized selective channel is considered, we 
obtain the same average timing metric. 
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gives rise to a carrier-dependent phase rotation at the fast Fourier 
transform (FFT) output, which can easily be compensated by 
the equalizers. Further, simulation results on the average timing 
metric are shown. As can be observed, the analytical results well 
correspond with the simulation results. Only at very low SIR, 
the analytical curve starts to deviate slightly from the simula- 
tion results: for very small SIR, the approximations used in the 
theoretical analysis are no longer valid. The narrowband inter- 
ference causes a decrease of the maximum of the timing metric 
when the SIR decreases, i.e., when the power of the interferer 
increases. As a result, the risk that the timing metric will not ex- 
ceed the threshold will increase such that the probability that the 
receiver will conclude that there is no training symbol present 
will increase. 

Fig. 6 shows results for the mean and variance of the timing 
metric at the optimum timing instant, as function of the SIR. For 
a large range of SIR, the theoretical results and the simulation 
results agree well. At large SIR, the mean of the timing metric 
reaches an asymptote, corresponding to the case where no inter- 
ference is present. This asymptote depends only on the SNR and 
is given by (1 + l/SNR)-^. At very low SIR, the mean of the 
timing metric becomes independent of the SNR. The theoretical 
results give rise to a maximum timing metric that is smaller than 
for the simulations when the SIR is small. Hence, the theoretical 
results can be considered as a lower bound on the performance 
for small SIR. The variance of the timing metric decreases for 
increasing SIR. This can easily be explained as at high SIR and 
SNR, the effect of the noise and interference is small, such that 
the symmetry of the pilot symbol is almost not affected, and the 
timing metric will be close to its average. 

Fig. 7 shows results for the mean and variance of the metric at 
the optimum timing point as function of the interference band- 
width. The results indicate that at given SIR, the mean is inde- 
pendent of the interference bandwidth. Further, the dependency 
of the variance on the interference bandwidth increases when 
SIR decreases. This is explained as at high SIR, the effect of in- 
terference diminishes and the variance mainly depends on noise. 
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At low SIR, the variance decreases when the interference band- 
width increases. 

Fig. 8 illustrates results for the mean and variance of 
the timing metric at the optimum timing point as func- 
tion of the normalized interference carrier frequency 
deviation F. The shape of the mean may be explained 
with the aid of Fig. 2. As long as the interfering signal 
is located in the flat region of the OFDM spectrum, i.e., 
ir-K(l + aO/(l + ao))(ro/T()| < ((l-ao)/(l + ao)), 
the mean is constant and does not depend on F. When 
the interfering signal is located outside OFDM band- 
width, i.e., (F + ((l-j-Q()/(l + ao))(To/T() > 1 or 
r - {{l + ai)/{l + ao)){To/Ti) < -1), the interfering 
signal does not affect the OFDM signal anymore: in this region 
the mean is constant and its value depends on the SNR value 
only. The periodicity of T( (see Fig. 3) does not appear in Fig. 8 
because Tj is much smaller than ^i. Therefore, the contribution 
of Ti on the mean and variance of the timing metric can be 
neglected. 
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'A' : SIR decreases as number of interferers increases 
' : SIR is constant for any number of interferers 
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Fig. 11. Mean and variance of timing metric outside training sequence at 
(SNR = 10 dB, NBW = 0.0244, T = 0, and = 1). 
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Fig. 9 shows the mean and the variance of the timing metric 
at the optimum synchronization point as function of the number 
of interfering signals, Nj, in two cases. In case "A," we consider 
that the SIR is fixed per interferer, so the total SIR decreases in- 
versely proportional to Nj . In case "B," we consider a fixed total 
SIR, i.e., SIR per interferer decreases linearly as Nj increases. 
As the mean and variance mainly depends on the total SIR, it 
follows that the mean and variance will be decreasing functions 
of the number Nj of interferers in case "A," whereas they are 
independent of Nj in case "B." 

Let us consider the case of two interferers (Nj = 2). We de- 
fine the SIR for the first interferer (SIRi) as P ■ SIR. Noting 
that (6) may be rewritten as SIR = l/Eti 1/SIRi, where 
SIRj = {{2a^/To)/{Ei/Ti)), the SIR for the second interferer 
(SIR2) is SIR/(1 - 1/(3). Fig. 10 shows the mean and the vari- 
ance of the timing metric at the optimum synchronization point 
as function of /3. The results indicate that the statistical prop- 
erties of the timing metric at the optimum timing point do not 
depend on SIR per interferer but only depend on the total SIR. 



Fig. 1 1 illustrates the results for the mean and variance of 
the timing metric at a timing instant outside the ti-aining symbol 
as fimction of SIR. At high SIR values, the interference signal 
does not affect the statistical properties of timing metric outside 
training symbol. Note that the simulation results are slightly dif- 
ferent from the analytical results. This is caused by correlations 
between the signal terms r* {{m+d)To) and r{{m+d+N/2)To) 
outside the training symbol, while the theoretical analysis as- 
sumed that they were independent. 

Fig. 12 displays the results for the probabiUty of missing a 
training sequence Pmiss at threshold A = 0.4.* The results show 
that Pmiss strongly depends on the NBI signal when the SIR is 
less than 15 dB. This is explained as the NBI signal destroys 
the correlation between the two halves of the received pilot 
sequence. 

On the other hand, Fig. 13 illustrates the analytical results of 
the probability of false detection Pfake as a function of SNR at 

■•The threshold level should be adapted according to the desired probability 
of missing a training sequence and probability of false detection. 
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different values of SIR. Note that at SIR = oo dB, the chi-square 
distribution is not a function of SNR, but only of the length of 
the training sequence. This explains why Pfaise is independent 
of the SNR at SIR = oo dB. Further, the values of Pfake are 
very small, approximately equal to zero. This can be explained 
as the length of the correlation (N = 1024) is large enough to 
cancel out the effect of noise, OFDM data, and the NBI signal. 
It is clear from the figure that the probability of a false detection 
is too small to be checked by computer simulations. 

V. Conclusion 

This paper evaluates the performance of the S&C symbol 
timing estimator [1] in the presence of narrowband interference 
for different interference signal parameters. The statistical prop- 
erties of the estimator are evaluated analytically and by means 
of simulation at the optimum synchronization point and outside 
the training sequence. Further, the probabilities of both missing 
and false detection of the training sequence are investigated. The 
agreement between the theoretical and simulation results proves 
the validity of our analysis. Generally, results indicate that the 
estimator works well for a wide range of SIRs. Also, at high SIR 
values (> 10 dB), noise is the dominant factor on the statistical 
properties of the estimator. Furthermore, at a given value of SIR, 
the interference bandwidth does not influence the mean of the 
metric while the dependency of the variance of the metric on the 
interferer bandwidth increases when SIR decreases. Moreover, 
the interference carrier frequency deviation does not affect the 
statistical properties of the estimator as long as the interferer is 
located in the flat region of the OFDM spectrum. The statistical 
properties of the estimator at optimum timing point do not de- 
pend on the SIR per interferer but only depend on the total SIR. 
At high values of SIR, the statistical properties of the estimator 
at a position outside the training symbol are essentially inde- 
pendent of the interference signal. Moreover, the probability of 
missing a training symbol is strongly affected by NBI as long as 
SIR is less than 15 dB. Finally, the probability of false detection 
approximately equals zero independent of the value of SNR and 
SIR, 



Appendix 
a, P, and 7 may be written as 
iv/2-1 

a= ^ {Re{r:(;jr„)r(3/„)}+j7m{r:(^„)r„(y„)}} 

(31) 

where Zm = {m + d)ro and = (m -f d -I- iNr/2)To. 

N/2-1 

^ Re{(rK^m) + «^(^n,))r(3/„.)} (32) 

N/2-1 

1= Re{(rJ(y„) + w;*(y„,))ru(s/n.)} 

m=0 

JV/2-1 

-|-2-R« ^ r*i{ym)-w{ym) 
N/2-1 

+ E {\w{ymt + \ri{ym)\'}. (33) 

fn=0 

After tedious computations, it follows that the mean and vari- 
ance of a, /3 and 7 can be written as 



N/2- 



04) 

(35) 



,.,=Re ^ ^(-'"')(0,f) 

m=0 ^ ^ 

+o/^\^"^'"^'\0,0) . (aC"-') (f>f ))* (36) 

-E'--"''(T-y) 

TT»,m'=0 ^ ' 

where B^"''"^'\x, y) is given by 

N, 

B^"^'"^\x,y) = Y, E {E[f-\m{{rr^+d+x)n-hTif 
■\gi{{m' + d^y)To-hTif. (39) 
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